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ABSTRACT 

A new method for the physical retrieval of rain rates from satellite microwave radiometers is 
presented and compared to two other rainfall climatologies derived from satellites. The method 
is part of a unified ocean parameter retrieval algorithm that is based on the fundamental 
principles of radiative transfer. The algorithm simultaneously finds near-surface wind speed W, 
columnar water vapor V, columnar cloud liquid water L, rain rate R, and effective radiating 
temperature Tu for the upwelling radiation. The performance of the algorithm in the absence of 
rain is discussed in Wentz [1997], and this paper focuses on the rain component of the algorithm. 

A particular strength of the unified algorithm is its ability to ‘orthogonalize’ the retrievals so 
that there is minimum cross-talk between the retrieved parameters. For example, comparisons of 
the retrieved water vapor with radiosonde observations show that there is very little correlation 
between the water vapor retrieval error and rain rate. For rain rates from 1 to 15 mm/h, the rms 
difference between the retrieved water vapor and the radiosonde value is 5 mm. A novel feature 
of the rain retrieval method is a beamfilling correction that is based upon the ratio of the retrieved 
liquid water absorption coefficients at 37 GHz and 19.35 GHz. This ratio decreases by about 
40% when heavy and light rain co-exist within the SSM/I footprint as compared to the case of 
uniform rain. This correction has the effect of increasing the rain rate when the spectral ratio of 
the absorption coefficients is small. Even with this beamfilling correction, tropical rainfall is still 
unrealistically low when the freezing level in the tropics (~ 5 km) is used to specify the rain layer 
thickness. We restore realism by reducing the assumed averaged tropical rain layer thickness to 
3 km, thereby accounting for the existence of warm rain processes in which the rain layer does 
not extend to the freezing level. 

Global rain rates are produced for the 1991 through 1994 period from observations taken by 
microwave radiometers (SSM/I) that are aboard two polar-orbiting satellites. We find that 
approximately 6% of the SSM/I observations detect measurable rain rates (R > 0.2 mm/h). Zonal 
averages of the rain rates show the peak at the intertropical convergence zone (ITCZ) is quite 
narrow in meridional extent and varies from about 7 mm/day in the winter to a maximum 1 1 
mm/day in the summer. Very low precipitation rates (< 0.3 mm/day) are observed in those areas 
of subsidence influenced by the large semipermanent anticylones. In general, these features are 
similar to those reported in previously published rain climatologies (Wilheit et al., 1991; Spencer, 
1993). However, significant differences do exists between our rain rates and those produced by 
Spencer (1993). These differences seem to be related to non-precipitating cloud water. 
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1. Introduction 

The accurate retrieval of oceanic rain rates from satellite microwave measurements has been 
an elusive goal since the concept was first proposed by Buettner (1963), and then demonstrated 
with satellite data by Wilheit et al. (1977). While the physical basis for these retrievals is sound, 
we believe that there are three significant assumptions inherent to these methods which make the 
measurement of tropical average rain rates to better than 10% problematic. These assumptions 
concern the specification of the following characteristics of a rain system: (1) the rain layer 
thickness (often assumed to extend from the surface to the freezing level); (2) the relative amount 
of cloud water versus rain water; and (3) the varying rain intensities across the radiometer 
footprint (which is commonly called the ‘beamfilling effect’). The observed brightness 
temperature (T B ) is strongly influenced by these three characteristics. Significant new 
information on these three issues will have to await the Tropical Rain Measurement Mission 
(TRMM, Simpson et al., 1988), scheduled to be launched in 1998. The combination of TRMM’s 
microwave radiometer and rain radar will help to quantify the above three processes. In the 
meantime, while the new rain retrieval method described herein does not solve these problems, it 
does attempt to explicitly address them in a physically realistic way. 

We present a unified, all-weather ocean algorithm that simultaneously finds the near-surface 
wind speed W, columnar water vapor V, columnar cloud liquid water L, rain rate R, and 
effective radiating temperature Tu for the upwelling radiation. This algorithm is a seamless 
integration of the Wentz (1997) no-rain algorithm and a newly developed rain algorithm. The 
algorithm is based on the fundamental principles of radiative transfer and explicitly shows the 
physical relationship between the inputs (T B ) and outputs (W, V, L, R, and Tu). The wind speed 
retrieval must be constrained to an a priori value for moderate to heavy rain, and Tu must be 
constrained by a statistical correlation for clear skies and light rain. The other retrieved 
parameters are unconstrained over the full range of weather conditions. Wentz (1997) discusses 
the algorithm’s performance in the absence of rain, and herein we focus on the rain component of 
the algorithm. 

A particular strength of the new method is its ability to ‘orthogonalize’ the retrievals so that 
there is minimum cross-talk between the retrieved parameters. With respect to estimating 
rainfall, it is important to remove the water vapor contribution to the observed brightness 
temperature. We will present results showing that the error in retrieved water vapor (as 
determined from radiosonde comparisons) is uncorrelated to the retrieved rain rate. Likewise, 
the influence of the radiating temperature Tu is separated from the liquid water signal by using 
the polarization information that is contained in the observations. Because the rain rates are 
retrieved only after all other significant influences on T B are quantified, the various retrievals can 
be analyzed for climate relationships between them, with high confidence that there is a 
minimum of algorithm cross-talk. 

Conceptually, the rain retrieval involves the following steps. The physics of radiative transfer 
shows that there is a direct and unique relationship between brightness temperature and the 
atmospheric transmittance Xl of liquid water. In view of this, the first step is to directly retrieve 
Tl along with the other directly observable parameters W, V, and Tu- In the context of rainfall, x L 
is related to the columnar water in the rain cloud, and Tu provides information on the height from 
which the radiation is emanating and whether radiative backscattering by large ice particles is 
occurring (Spencer, 1986). The retrieval of x L is done by solving a set of simultaneous brightness 
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temperature equations. A basic premise in this retrieval is that the polarization signature of the 
T b ’s allows for the separation of the x L signal from the Tu signal. The T B model is formulated 
such that the Tu parameter includes both radiative scattering effects and air temperature 
variability. In the next step, the spectral signature of the retrieved Tl at 19 and at 37 GHz is used 
to estimate the beamfilling effect. A beamfilling correction is applied, and the mean atmospheric 
attenuation A L for liquid water over the footprint is found. Mie scattering theory and an assumed 
relationship between cloud water and rain water are then used to convert the A L to a columnar 
rain rate, which is defined as the vertically averaged rain rate times column height. Finally, the 
columnar rain rate is converted to a vertically averaged rain rate by dividing by an assumed rain 
column height that is a function of a sea surface temperature climatology. The final assumption 
is that the surface rain rate equals the vertically averaged rain rate. In this way, we explicitly 
handle the three rain cloud characteristics listed above. 

The algorithm is developed and tested using the observations taken by the Special Sensor 
Microwave Imager (SSM/I, Hollinger et al., 1987). The SSM/I is a scanning radiometer that 
operates at four frequencies: 19.35, 22.235, 37, and 85.5 GHz. It is flown by the Defense 

Meteorological Satellite Program (DMSP) on two operational polar orbiting platforms. The first 
in the series of 7 SSM/I's was launched in June 1987. The results herein are based on SSM/I 
observations for the four year period from 1991 through 1994. Observations from two satellites, 
F10 and FI 1 , are used. The F10 observations cover the entire 4 year period, while the Fll 
observations begin in January 1992. The SSM/I time series will probably continue through about 
the year 2000, at which time it will be replaced by a combined imager/sounder called SSM/IS. 
Thus there is the opportunity to obtain a 13-year global time series of geophysical products, 
which can then be further extended with the SSM/IS observations. 

The algorithm described herein will be used to produce the NASA Pathfinder Data Set for 
SMMR and SSM/I. This data set will be a 20-year time series of geophysical parameters, which 
will be broadly distributed to the research community. 


2. Description of the SSM/I 

The SSM/I’s are carried aboard the DMSP polar orbiters. The orbit for these spacecraft is 
near-circular, sun-synchronous, and near-polar, with an inclination of 98.8°. The nominal 
altitude is 860 ± 25 km, and the orbital period is about 102 minutes. The variation in altitude is 
due to both the eccentricity of the orbit and the oblateness of the Earth. To date, five SSM/I’s 
have been launched aboard the F08, F10, Fll, F12, and F13 spacecrafts. Unfortunately, the F12 
SSM/I failed. Table 1 gives the launch dates and ascending node times for the remaining four 
SSM/I’s. The ascending node times for F10 and FI 1 are slowly drifting later into the day. The 
data from the F08 SSM/I essentially ends in December 1991. Data from the other three SSM/I’s 
are continuing as of December 1996. 

The SSM/I sensor consists of 7 separate total-power radiometers sharing a common feedhom. 
These 7 radiometers take dual-polarization measurements at 19.35, 37.0, and 85.5 GHz, and just 
a vertical-polarization measurement at 22.235 GHz. The SSM/I uses an offset parabolic reflector 
of dimensions 61 by 66 cm to collect the microwave radiation. The reflector focuses the 
radiation into the corrugated, broad-band, 7-port feedhom. The reflector and feedhorn spin as a 
unit about the nadir axis. The rotation period is 1.9 s. A cold-space reflector and a hot reference 
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load are attached to the spin axis and do not rotate. The rotating feedhom observes the fixed cold 
reflector and hot load once each scan. In this way, calibration observations are taken every scan. 


Table 1 . SSM/I Launch Dates and Ascending Node Times 


Satellite 

Launch Date 

Ascending Node Time 

F08 

June 1987 

6:15 AM 

F10 

Dec. 1990 

8:10 PM to 10:15 PM 

FI 1 

Nov. 1991 

5:10 PM to 6:30 PM 

F13 

March 1995 

5:45 PM 


Earth observations are taken during a 102.4° segment of the rotation. The 102.4° arc is 
centered on the spacecraft subtrack and corresponds to a 1400-km wide swath on the Earth's 
surface. The 1400-km swath and the orbit inclination of 98.8° provide complete coverage of the 
Earth in two to three days, except for two small circular sectors of 2.4° centered on the North and 
South poles. The nadir angle for the Earth-viewing reflector is 45°, which results in an Earth 
incidence angle of 53.4° ± 0.25°. The lower frequency channels (19, 22, and 37 GHz) are 
sampled so that the pixel spacing is 25 km, and the 85 GHz channels are sampled at a 12.5 km 
pixel spacing. More information on the SSM/I can be found in Hollinger et al. (1987) and Wentz 
(1991). 

3. The No-Rain Algorithm 

We begin by reviewing the no-rain algorithm described by Wentz (1997). Then Section 4 
shows how this algorithm is extended to include rain observations. In the absence of rain, there 
is a relatively simple and unique relationship between the ocean brightness temperature (T B ) 
measured by SSM/I and W, V, and L. As a consequence of this simple relationship, these 
parameters can be retrieved to a high degree of accuracy. The retrieval of (W,V,L) is 

accomplished by varying their values until the T B model function matches the SSM/I 
observations. After a precision calibration to in situ observations, the rms retrieval accuracies for 
W, V, and L are 0.9 m/s, 1.2 mm, and 0.025 mm, respectively (Wentz, 1997). We now give 
some details on the no-rain algorithm so that one can then see how it is extended to include rain. 

The first step in the retrieval process is to average the SSM/I observations to a common 
spatial resolution. The half-power beamwidths of the SSM/I footprints on the Earth are 56, 44, 
and 32 km for the 19, 22, and 37 GHz channels, respectively, and the centers of these footprints 
are coincident. Hence, a 37 GHz observation sees only 33% of the area sampled by the 19 GHz 
channel. In order to obtain accurate retrievals, it is necessary that all channels see the same ocean 
area. This is accomplished by averaging the 22 and 37 GHz observations down to the lower 
resolution of the 19 GHz channel using the following equation: 
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( 1 ) 


TaIJ ~ ^ S W ij^Aij 

i=I-lj=J-l 

in which T A ij is antenna temperature (either 22 or 37 GHz) at the original resolution and the two 
subscripts denote the along-track scan number and the across-track cell position. A set of 
weights Wjj (one for 22 GHz and another for 37 GHz) is found such that the effective antenna 
pattern of the averaged T A matches the 19-GHz pattern. The weights depend on the across-track 
cell position because the relative location of the cells is different at the swath edge as compared 
to the swath center. The distance between adjacent scans and adjacent cells is approximately 
25 km, and we find that it is sufficient to include only the immediately adjacent cells when doing 
the average. Hereafter, we drop the overbar on T A , and it is understood that all observations are 
at a common spatial resolution. Thus, the retrievals are done at a pixel spacing of 25 km, but 
their actual spatial resolution is about 56 km. Note that for the rain algorithm, the rain rates are 
computed at the spatial resolution of the 37 GHz footprint, which is 32 km (see Section 6). 

The antenna temperatures are then converted to brightness temperatures using the linear 
transformation given by Wentz (1991). The T A -to-T B conversion accounts for the antenna 
spillover and the cross-polarization leakage. The no-rain algorithm then retrieves three 
geophysical parameters over the ocean: near-surface wind speed W (m/s), columnar water vapor 
V (mm), and columnar cloud liquid water L (mm). The retrieval is done by solving the following 
three model function equations: 

T A22 v = G V v[F 2 2 v(W,V,L) + 0.5A(x)x 2 2 AT BV ] 

+ GvhF 22 h(W,V,L) + GvoTbc (2a) 

TB37V = F 37 v(W,V,L) + 0.9A(x) x] 7 AT BV (2b) 

T B 37h = F 37 h(W,V,L) (2c) 

The SSM/I brightness temperatures at 37 GHz are denoted by T B37 v and T B37 h for vertical and 
horizontal polarization, respectively. At 22 GHz, there is only vertical polarization, and the 
model function equation is in terms of the antenna temperature T A 22 v rather than brightness 
temperature. The term T B c is the cosmic background brightness temperature equaling 2.7 K. 
The G factors are the antenna pattern coefficients that account for antenna spillover and cross- 
polarization leakage. The values for the G coefficients are given by Wentz (1991). The function 
F is the Tb model function for an isotropic sea surface (i.e., no wind direction effects). In 
addition to the explicit arguments W, V, and L, F also depends on sea surface temperature T s . A 
sea-surface temperature climatology (Shea et al, 1990) is used to specify Ts- 

The terms containing AT B v account for the dependence of T b on wind direction. AT B v is the 
difference between the 19 GHz SSM/I v-pol brightness temperature T B i 9 v minus the isotropic 
model function and hence is a measure of wind direction variations. 

AT B v = T B i9v — F|9v(W,V,L) (3) 

The wind direction signal is diminished by the square of the transmittance, x 2 . At low winds 
for which the wind direction signal is weak, the modeling noise in the retrieval equations is 
reduced by introducing the term A(x): 
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A(x) = 0 

x < 0 

(4a) 

A(x) = 3x 2 - 2x 3 

0 < x < 1 

(4b) 

A(x) = 1 

x > 1 

(4c) 

x = (W - 3)/5 


(5) 


A(x) is a weighting function that smoothly goes from 0 to 1 as its argument x goes from 0 to 1 . 
The coefficients 0.5 and 0.9 in equations (2a) and(2b) are chosen so as to minimize the rms 
retrieval error. 

Equations (2) represent three equations in three unknowns: W, V, and L. These three 

equations are solved by assuming the equations are stepwise linear in terms of (W,V,L). The 
analytic derivatives 3 Tb/3W, 9Tb/ 9V, and dT^/dL are computed at a first guess value for 
(W,V,L). The set of equations is then treated as a linear system with slopes equal to the first 
guess analytic derivatives, and this set of equations is solved in the usual way using Cramer's 
rule. The solution gives a new estimate for (W,V,L). The analytic derivatives are then 
recomputed at the new solution point, and the equations are again solved. This procedure is 
continued until the difference between the observation and the model function is less than 0. 1 K 
for each channel. Typically 5 iterations are required to reach the 0.1 K convergence level. 
Although this iterative procedure requires a first guess, it should be emphasized the final solution 
is independent of the first guess because equations (3) and (4) are deterministic: there are 4 
independent equations and 4 unknowns. 


4. Extending the Algorithm to Include Rain 

To create an all-weather algorithm, the no-rain algorithm is extended in the following ways: 

1 . The cloud water parameter L is replaced by the total transmittance of cloud and rain 
water at 37 GHz, Tu 7 - 

2. An additional parameter is added to the retrieval, the cloud and rain transmittance at 19 
GHz, Tli 9 - 

3. When rain occurs, the wind speed retrieval is constrained to an a priori value. 

4. When rain occurs, the effective air temperature Tu for the upwelling radiation becomes 
a retrieved parameter. For the no-rain algorithm, Tu is specified as a function of 
retrieved water vapor and climatological sea-surface temperature (Wentz, 1997). 

The algorithm extensions are done such that the all-weather algorithm is identical to the no-rain 
algorithm when there is no rain, and then smoothly transforms to a rain algorithm over the range 
from very light to moderate rain. Above about 2 mm/h the algorithm no longer retrieves wind 
speed. The first three modifications are discussed in this section, and the fourth modification 
involving the effective air temperature is discussed in Section 6. 

The first modification is to replace the argument L by the 37-GHz liquid water transmittance 
Tl 37 for both cloud and rain water. In the no-rain algorithm, the cloud liquid water L enters the 
equation only through the transmittance. This dependence can be expressed as 

x = exp [-sec 0 (A 0 + A v + A L )] (6) 

A L = <Xl [1 - 0.026(T l - 283)] L (7) 
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where 0 is the incidence angle, A 0 and A v are the atmospheric oxygen and water vapor 
absorptions, and A L is the cloud liquid water absorption. Wentz (1997) gives expressions for A 0 
and A v as functions of the effective air temperature and columnar water vapor. In the absence of 
rain, the radiative transfer through the cloud droplets, which are much smaller than the radiation 
wavelength, is governed by Rayleigh scattering, and the absorption is proportional to the 
columnar liquid water content L (mm) of the cloud (Goldstein, 1951), as indicated by (7). The 
Rayleigh absorption coefficients (Xl are 0.059, 0.078, and 0.208 for 19, 22, and 37 GHz, 
respectively, when L is in units of millimeters. The absorption has a small dependence on the 
temperature T L of the cloud water, which is assumed to be the mean temperature between the 
surface and the freezing level. 

The liquid water dependence in (6) can be equivalently expressed in terms of the liquid water 
transmittance Tlf: 

t f = t L F exp[-sec 0 (Ao + A v )] (8) 

t lf = exp(-A L Fsec 0) (9) 

The subscript F has been introduced to denote the dependence on frequency. Often it is more 
convenient to work in terms of the absorption than the transmittance. Given the retrieved value 
for Tlf. the absorption (before doing the beamfilling correction) is given by 

Alf = -cos 0 ln(TLF) (10) 

When rain is present, the relationship between Tlf and liquid water content is more complex, 
as discussed in Section 8, and the simple Rayleigh expression is not valid. However, by 
parameterizing the T B model in terms of Tlf rather than L, we defer the problem of relating Tlf 
and liquid water content. In other words, we are dividing the rain retrieval problem into two 
steps. This first step involves separating the liquid water signal, expressed in terms of Tlf. from 
the signal of the other parameters. Since T B is nearly proportional to Tlf 2 (see equation 29), the 
beamfilling problem is not an issue when retrieving Tlf- In the second step, a rain rate is inferred 
from Tlf- It is only in the second step that one is required to make assumptions regarding the 
beamfilling, the cloud/rain partitioning, and the rain column height. 

The second modification is to introduce T L i 9 as an additional parameter to be retrieved. For 
the no-rain algorithm, Rayleigh scattering gave a fixed relationship between the 19, 22, and 37 
GHz transmittances, and hence it was not necessary to separately retrieve Tli 9 - However when 
rain is present there is no fixed spectral relationship between the transmittances. Accordingly, 
we directly retrieve Tli 9 by introducing a fourth equation in to the retrieval process. Thus, Tli 9 
and Tl 37 are retrieved products. The 22-GHz transmittance t L 22 is found from an interpolation 
between Tli 9 and Tl 37 , as is discussed below. 

The third modification is to eliminate wind speed as a retrieved product when there is 
significant rain. The decrease in the atmospheric transmittance obscures the surface and 
degrades the ability to retrieve the wind speed. Furthermore, the T B modeling error is larger for 
raining observations due to errors in specifying the effective air temperature, as is discussed in 
Section 6. For moderate to heavy rain it is best to constrain the wind parameter to some specified 
a priori value. To do this, we use the SSM/I wind retrievals in adjacent, no-rain areas to specify 
W. If no such wind retrievals are available, we use a monthly, 1° latitude by 1° longitude wind 
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climatology to specify W. This climatology is produced from 7 years of SSM/I observations. 


The wind direction effect (i.e., AT B v terms in (2)) only weakly affects the retrievals and is set to 
zero for raining observations. 

These three modifications result in the following set of retrieval equations: 

T Bl9 v = F 19 v(W,V,t l , 9 ) + 0.9A(x)[ l-A(y)] tf 9 AT BV (11a) 

T A 22 v = Gw { F 22 v(W ,V. ,Ti_ 22 ) + 0.5A(x)[ 1 — A(y)] x 22 AT B v } 

+ Gvh F 22 h(W,V,Tl 22 ) + GvoT B C (11b) 

T B3 7 V = F 37 v(W,V,t L 37 ) + 0.9A(x)[l— A(y)] t’ 7 AT bv (11c) 

Tb37H = F 37 h(W,V,T L 37) - My) [F37H(Wo,V,TL37) - T B37H ] ( 1 1 d) 

The constraint on wind is accomplished by the A(y) term, which is the weighting function given 
by (4) above. In this case the argument y is 

y = (A L37 - 0.04 )/0. 1 1 (12) 

Thus A(y) smoothly goes from 0 to 1 as its argument Al 3 7 goes from 0.04 to 0.15. The lower 


bound for Al 3 7 , which represents the onset of rain, was found from an investigation of 38 
northeast Pacific extratropical cyclones (Wentz, 1990). This study indicated that when A L37 
exceeds 0.04, drizzle or light rain is likely. The upper limit represents a rain rate of about 
2 mm/h, depending on the rain column height and other factors. When A(y) = 0, the above 
equations have the same form as the no-rain equations (2). When A(y) = 1 , the wind direction 
term is eliminated, and equation (lid) simply becomes 

F 37H (W,V,t L37 ) = F 37H (W 0 ,V,Tu7) (13) 

where Wo is the a priori wind speed. Equation (13) forces the wind retrieval to equal Wo. Thus 
Al 37 plays the role of a blending parameter for joining the no-rain algorithm with the raining 
algorithm. As A L37 goes from 0.04 to 0.15, the algorithm smoothly transforms from the no-rain 
algorithm to a rain algorithm. 

In the no-rain algorithm, Rayleigh scattering gives the following relationship between the 22 
and 37 GHz transmittances: 


_ _ „0 375 

^L22 — ^L37 


(14) 


When rain is present, the Rayleigh relationship is no longer valid, and we use the following 
interpolation based on the retrieved values of the 19 and 37 GHz transmittances: 


T 


L22 


0.837 

T L19 


0.163 

l L37 


(15) 


This interpolation represents a linear interpolation of A L versus frequency f. According to Mie 
scattering, the spectral variation of A L is approximately f 1 7 for high absorption values. However, 
the beamfilling effect discussed in Section 7 tends to flatten out the spectral variation of A L , and 
a linear variation is found to be a good approximation for moderate to high rain rates. In order to 
blend the no-rain algorithm with the raining algorithm, the following expression is used to 
specify T i ■>? for all values of A L37 
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t l22 = [(1 - A(z)]t° l ™ + A(z)<ff ( 1 6) 

Z = 10Al37 • (17) 

Thus for small A L 37 , Tl 22 is specified via Rayleigh absorption, and for large A L37 , x L22 is specified 
via the 19-37 GHz interpolation. 

The remaining equation to be modified is equation (3) for the 19V T B wind direction signal 
AT B v- In this equation, the L dependence is specified in terms of the 19 GHz liquid water 
transmittance Xu 9 - Since this wind direction term is set to zero for rain rates above about 2 
mm/h, we can use the Rayleigh relationship to specify Tli 9 as a function of Tu 7 - Doing this, we 
obtain the following expression for AT B y: 

AT B v — T B 19V — Fl9v(W,V,T°37 6 ) (18) 

where the Rayleigh transmittance at 19 GHz. 

With the above modifications, there are now four equations in four unknowns (W, V, Tli 9 . and 
'lu.i)- This set of equations is solved in the same way as was discussed above for the no-rain 
algorithm. The equations are assumed to be stepwise linear in terms of (W,V,Tli 9 Al 37 X an d the 
equations are solved in an iterative manner. Typically about 5 iterations are needed to reach the 
0.1 K convergence level. 

In summary, when it is not raining and A 137 is small, the extended algorithm is identical to the 
no-rain algorithm described by Wentz (1997). Then as A L37 increases, the algorithm smoothly 
transforms to a rain algorithm. The wind estimate is constrained to the a priori value, and the 
fixed Rayleigh relationships among the 19, 22, and 37 GHz transmittances are relaxed, but the 
same T B model function is used throughout. In this way, we obtain a unified, all-weather ocean 
algorithm for SSM/I. 


5. Retrieval of Water Vapor in Rain 

Figure 1 shows the difference between the SSM/I retrieved water vapor and the value obtained 
from collocated radiosonde observations (RAOB). The difference is plotted versus rain rate. 
The quality control of the radiosonde data and the collocation with the SSM/I are discussed in 
Wentz (1997). There are a total of 35,108 SSM/I overpasses of radiosonde sites. For these 
overpasses, a total of 81 ,922 rain observations are found within a 1 12 km radius of the site. The 
solid curve shows the mean difference and the dashed curves show the ± one standard deviation 
of the difference. The rain rate is computed from the SSM/I observations, as described in this 
paper. The statistics are computed by first binning the observations into 0.5 mm/h rain rate bins. 
For rain rates between 1 and 15 mm/h, the typical rms difference between the SSM/I and 
radiosonde vapor is 5 mm. In comparison, the rms difference for the no-rain observations is 
3.8 mm. The error analysis in Wentz [1997] indicates that the spatial and temporal sampling 
mismatch between the SSM/I 56-km footprint and the radiosonde point observation contributes 
about 3.7 mm to the total rms difference. Thus nearly all of the rms difference for the no rain 
observations is due to the spatial-temporal mismatch. For the rain observations, about half of the 
rms difference is due to the spatial-temporal mismatch. The conclusion we draw from Figure 1 


8 



is that the retrieval algorithm does provide reasonably good estimates of water in the presence of 
rain and there is little cross-talk between the vapor retrieval and the rain retrieval. 



Fig. 1 . A comparison of SSM/I and radiosonde columnar water vapor for rainy 
observations. The solid line is the mean difference, and the dashed lines show 
the one standard deviation envelope for the differences. 


When we first compared the SSM/I and radiosonde water vapor, we found that for large A U9 , 
the SSM/I retrievals were biased low (high) relative to the radiosonde values for high (low) 
values of sea-surface temperature T s . In view of this, we apply the following ad hoc correction 
to the retrieved values of V in order to correct the observed systematic error. 

Vcorr = v ret + 12 (Ts - 283) (Alw - 0.05)' 5 (19) 

This correction has been applied to Figure 1. We believe that this systematic error is due to 
radiative scattering and the degradation in the assumed relationship between the effective air 
temperature and water vapor, as is discussed in the next section. For rain rates above about 
10 mm/h, the size of this correction is about 10 mm. 

We find that when A L | 9 exceeds about 0.3 (which corresponds to R - 15 mm/h, depending on 
rain column height), the atmosphere is too opaque and/or scattering is too strong to obtain a 
useful estimate of V. The procedure discussed in Section 6 for obtaining tL when radiative 
scattering is significant requires that V be specified. Thus, for A L i 9 > 0.3, we use an a priori 
value for V based on the SSM/I vapor retrievals in adjacent, no-rain and light-rain areas. If no 
adjacent retrievals are available, we use a monthly, 1° latitude by 1° longitude vapor climatology 
to specify V. This climatology is produced from 7 years of SSM/I observations. 
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6. Effective Air Temperature and Radiative Scattering 

Implicit in the T B model function F(W,V,t l ) is an assumed relationship for the effective air 
temperature Tu versus the retrieved columnar water vapor V and a climatological sea-surface 
temperature Ts- The effective air temperatures, which we also call the effective radiating 
temperatures, for the upwelling and downwelling atmospheric radiation are defined as 

Tu = T B u/(1-t) (20a) 

T D = T BD /(l-t) (20b) 

where T B u and T B d are the upwelling and downwelling atmospheric brightness temperatures and 
x is the atmospheric transmittance. An analysis of 42,195 radiosonde observations shows that in 
the absence of rain Tu and T D are well correlated with the V and T s (Wentz, 1997): 

Tu = ^(V,T S ) (21) 

T d = Tu-c 6 -c 7 V (22) 

where the function H^V/Ts) and the coefficients C6 and c 7 are given by Wentz (1997). At 19 and 
37 GHz, Tu and T D are very similar, with Tu being about 2 K colder. 

When rain is present, the accuracy of 4'(V,T S ). degrades due to two factors. First the process 
of precipitation and associated convection alters the air temperature. Second, radiative scattering 
by large rain drops and ice particles reduces T B u and T B d- Since by definition Tu and Tp are 
proportional to T B u and T B d, they also decrease when scattering occurs. Fortunately, the 
radiosonde comparisons in the previous section show that the retrieval error in V due to the 
degradation of 'F(V,T S ) is not large and can be partially corrected by the ad hoc adjustment (19). 
However, the error in Xu 9 and x L 37 can be large if scattering effects are ignored. 

The effect of radiative scattering on the T B model is only through the atmospheric parameters 
x L , T bu , and T BD - In view of definitions (20), the T BU and T BD dependence can be equivalently 
expressed in terms of Tu and T D . The rain algorithm accounts for scattering (and rain-induced 
variations in air temperature) by letting Tu become a retrieved parameter rather than being 
specified by ^(V.Ts). Since V and W have been already found by the procedure described in 
Section 4, the retrieval problem is reduced to the following two equations in two unknowns (i.e., 
Tu and x L ). 

T bv = F v (Tu,x l ) (23a) 

T bh = Fh(Tu.Xl) (23b) 

There is a separate pair of equations for 19 and 37 GHz, and we have temporarily dropped the 
subscript denoting frequency. Two assumptions are required to solve these equations. First, Tu 
and To are assumed to be closely correlated so that To can be specified as a function of Tu, as 
given by (22) above. The second assumption is that Tu has the same value for vertical and 
horizontal polarization. In the absence of scattering, Tu is completely independent of 
polarization. And for moderate to heavy rain, for which scattering is significant, SSM/I 
observations show that the saturation values for the v-pol and h-pol T B observations are nearly 
the same (Spencer et al., 1989). Thus the assumption of polarization independence seems 
reasonable. 
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Section 7 gives an example of how equations (23) are inverted to yield Xl and T(j assuming a 
simplified T B model. We denote the solution for Tu for the complete T B model by the function 

Tu = ^(T bv , T bh ) (24) 

We use this expression to specify the effective air temperature for the case of moderate to high 
rain rates. In order to blend this estimate of T y with the no-rain estimate given by 4'(V,T S ), we 
use the following expression: 

Tu = [ 1 - A(u)] y(V,T s ) + A(u) ^(T bv , T bh ) (25) 

u = (0.7 - t)/0.2 (26) 

and x is the total transmittance x of liquid water, water vapor, and oxygen. Equations (23a) and 
(25) are then combined to retrieve x L . When x > 0.7, then Tu = H'fVJ's), and the retrieved x L is 
identical to that found by the no-rain algorithm described in Section 3. When the x < 0.5, the 
retrieved x L and Tu are based solely on the magnitude and polarization signature of the dual- 
polarization observations (either 19 or 37 GHz). In this way, the all-weather algorithm is 
identical to the no-rain algorithm for x > 0.7. Then as x decreases, the equations smoothly 
transform to a rain algorithm that characterizes the emission and scattering by rain via the Xl and 
Tu retrievals, which are based just on the signature of two, dual-polarization observations that 
have been corrected for water vapor and wind speed effects. Note that £,(T B v. T BB ) is inversely 
proportional to 1 - x, and large retrieval errors occur when x ~ 1. Thus we only use this method 
when x < 0.7. Figure 2 below shows how the accuracy of £(T BV , T BH ) degrades for x > 0.7. 

We define the depression ATu in the effective temperature as 

ATu = ^(T bv , T bh ) - ^(V,T S ) (27) 

which is a measure of the decrease in brightness temperature due to radiative scattering. In 
addition to radiative scattering, ATu is also a measure of the decrease in the air temperature due 
to most of the radiation coming from the cold cloud tops. Figure 2 shows ATu plotted versus x 
for the time period from July through September 1992, as derived from the F10 observations. 
For this time period, the SSM/I retrieval algorithm finds 7,859,295 rain-influenced footprints 
over the world’s oceans. The curves are generated by binning these observations into Ax = 0.01 
bins. The solid curves show the mean value for each bin, and the dashed curves show the ± one 
standard deviation for each bin. The depression in the effective temperature for the 37 GHz 
observations, which are most affected by scattering, is about twice that of the 19 GHz 
observations. For low values of x, ATu is about —10 K and —20 K for 19 and 37 GHz, 
respectively. For x > 0.7, the retrieved value of Tu given by ^(T BV ,T BH ) for 19 GHz becomes 
noisy and unreliable. 

The retrievals W, V, x L and Tu are all done at the common spatial resolution of the 19 GHz 
channels, which is about 56 km, as dictated by equation (1). For the rain rate retrievals, we want 
as much spatial resolution as possible. In order to obtain a rain rate at the resolution of the 37 
GHz footprint, we make the assumption that W and V are uniform over the 19 GHz footprint. 
The above equations are then used to find XL37 and Tu37 given the 37-GHz T B s at their original 
resolution of 32 km. In Figure 2, the spatial resolution for the ATu values is 56 km for 19 GHz 
and 32 km for 37 GHz. 
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Fig. 2. The decrease in the effective air temperature due to 
radiative scattering and cold cloud top temperatures. 


The calculation of beamfilling correction factor discussed in the next section requires both 
Tli 9 and Xl 37 - Thus to compute a rain rate at the 32-km resolution, a value of Xli 9 at the spatial 
resolution of 32-km is required. We use the following expression to specify a high-resolution 

XL19.HI 

ln(TLi9,Hl) = ln(XLI9) [ln(XL37.Hl)/l n (XL37)] (28) 

where the subscript HI denotes the higher spatial resolution. Note that (28) is in terms of the 
absorption, i.e., ln(x L ). The assumption behind equation (28) is that, even though the cloud and 
rain water may significantly vary over the footprint, the ratio of 19-to-37 GHz liquid water 
absorption is relatively constant. If this is true, then the observed spatial variation in x U 7 can be 
used as a scaling factor to compute Xli 9 at the higher spatial resolution. In reality, this spectral 
ratio will have some inter-footprint variability, and (28) will introduce some error into the rain 
retrieval. However, this error will tend to be unbiased and will tend to zero when doing spatial 
and temporal averages of rain rate. 
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7. The Beamfilling Effect 

The terminology ‘beamfilling effect’ refers to the non-linear relationship between rain rate R 
and brightness temperature T B that occurs when averaging over the radiometer footprint. To 
illustrate this effect, we use a relatively simple model for the brightness temperature: 

T b = T e (t)(1-t 2 p) (29) 

where T E (x) is the effective temperature of the combined ocean and atmosphere system, x is the 
total transmittance through the atmosphere, and p is the reflectivity of the sea surface. The 
effective temperature is a function of x. For x = 1, T E equals the sea surface temperature T s , and 
for x = 0, T e equals the effective temperature Tu of the upwelling atmospheric radiation. T E 
smoothly varies from T s to Tu as x goes from 1 to 0. 

Despite its simplicity, equation (29) is actually a very good approximation to the complete T B 
model function described in Section 3 and provides considerable insight into the retrieval 
algorithm discussed in the previous sections. As discussed in Section 6, to separate x from the 
T e signal requires dual-polarization observations, either at 19 or 37 GHz. Looking at the simple 
T b model above, we see that the T e term is easily eliminated, and the transmittance is given by 

t 2 _ T bv~ T bh — (30) 

Ph^BV — Pv^BH 

Values for the wind speed and water vapor come from the procedure discussed in Section 4. The 
wind speed and climatology sea-surface temperature are used to specify pv and p H . The oxygen 
and water vapor components of x 2 are factored out using equation (8), thereby obtaining just the 
liquid water transmittance Xl 2 . This method of using polarization information to separate the x 
signal from the T E signal is similar to that proposed by Petty (1994). 

We see that for rain retrievals the basic observable is x L 2 , which represents an average value 
over the SSM/I footprint. Using equation (9), the footprint average is given by the following 
integral 

x 2 L =|e" 2A ' sec0 P(A , )dA / (31) 


where the integral is over the probability density function P(A') for the cloud and rain water 
absorption A' within the footprint. The desired retrieval is the mean absorption over the 
footprint, which is given by 

A L = J A'P(A') dA' (32) 

Without any correction for the spatial averaging, the estimate of the absorption would be 


A = ~ lntL 
L sec0 


(33) 


Any variation of the absorption within the footprint will result in the estimate Al being lower 
than the true mean value A L . This systematic underestimation of the absorption is called the 
beamfilling effect. 


If P(A') is an exponential distribution of the form 
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34) 


P(A') 



then the integrals in (3 1) and (32) are easily evaluated, and one obtains the following relationship 
between Al and A L : 


g2A L p 2 sec0 _ j 

2p 2 sec0 


(35) 


P = 


aa l 

A l 


(36) 


AA l is the rms variation of A' within the footprint and P is the normalized rms variation of A'. 
For an exponential distribution p equals unity. An expansion of statistical moments to the 
second order in A L p 2 shows that (35) is correct for any distribution P(A') having a mean A L and 
an rms variation of AA L . Thus if P is known, (35) can be used to correct the beamfilling effect. 
P is a normalized quantity that is related to the variability of liquid water in the footprint and 
hence is essentially independent of frequency. 

In order to find p, we note that the ratio Al 37 /Ali 9 will be less than 
that predicted by Mie scattering theory when the beamfilling effect 
is significant. Thus comparing A^/Aug to the Mie absorption 
ratio provides the means to determine the beamfilling effect. The 
relationship between the Mie ratio and Al 37 /Ali 9 is given by the 
37-to-19 GHz frequency ratio of equation (35): 


Al37 _ 
Al19 


e 2A L37 p 2 sece _ j 
g 2 ^L19p sec ® J 


(37) 


where the left-side is the Mie ratio given by equation (40) below. This ratio varies from 3.5 for 
light absorption values to 2.8 for heavy absorption. If the beamfilling effect is not significant, 
then Al 37 /Ali 9 will equal Aijn/Auy. (Because of retrieval error, Al37/Ali 9 is sometimes slightly 
greater than A U 7 /A L i 9 .) Thus when A L 3 7 /A L i 9 > A^/Alu, no beamfilling correction is done 
(i.e., p is set to 0). Otherwise (37) is inverted to find p as a function of the two ratios: Al37/Ali 9 
and Al 37 /A L i 9 . The inversion is readily done because the expression on the right-side of (37) 
increases monotonically with p. It equals Al 37 /Ali 9 when p = 0. Thus, the beamfilling 
correction consists of finding a value for P that when substituted into (37) yields the absorption 
ratio given by the Mie theory. Given p, the mean absorption for the footprint is then found (35). 
Since the Mie ratio depends on Al 37 , an iterative procedure is required to find p. However, the 
iterative process quickly converges to a unique solution for p because A L 37 /Ali 9 only weakly 
depends on Au 7 - 

The magnitude of the beamfilling correction is characterized in terms of the ratios A L j 9 /A U 9 
and Al 37 /Al 37 > which are called the beamfilling correction factors (BCF). When A L 3 7 /A U9 is 
significantly less than A L 37 /A L i 9 , large values for p and BCF are found. For example, when 
A U7 /A L i 9 = 2, the BCF is 1 .4 and 2.0 for 19 and 37 GHz respectively. For even smaller values of 
A L 37 /Ali 9 , the BCF increases exponentially, and we must impose the following limits. The 
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maximum values of 3.4 and 6.4 are used for the 19 and 37 GHz BCF, respectively, which 
corresponds to the exponent 2Au7p 2 sec0 in (35) reaching a value of 3.0. If the BCF exceeds the 
maximum, it is reset to the maximum. Another overall limit is placed on A L | 9 and A L3 7 - Neither 
values is allowed to exceed 1.2. These limits correspond to observations of heavy rain for which 
the 37 GHz and, to a lessor degree, the 19 GHz brightness temperatures have reached saturated 
levels. The effect of these limits is to place an upper bound on the retrieved rain rate. For the 
extreme case of A U9 reaching a value of 1.2, the retrieved rain rate will be about 25 ram/h (75 
mm/h) for a rain column height of 3 km (1 km). We consider the 25 mm/h limit as an extreme 
upper bound on the algorithm’s ability to retrieve rain. For such high rain rates, both the 19 and 
37 GHz observations have saturated, and the retrieval error can be very large. A more 
conservative approach would be to lower the limit to 15 or 20 mm/h (for a 3 km column height). 

Figure 3 shows the 37-GHz absorption plotted versus the 19-GHz absorption for the July 
through September 1992 period discussed in Section 6. The bottom curve in Figure 3 shows the 
retrieved absorptions Aurz versus A L | 9 before the beamfilling correction. The middle curve 
shows the absorptions A L37 versus A L!9 after the beamfilling correction, and the top curve shows 
the theoretical curve computed from Mie scattering computations. The curves are generated by 
binning the 7,859,295 observations into absorption bins having a width of 0.005. The solid 
curves show the mean value for each bin, and the dashed curves show the ±1 standard deviation 
for each bin. The Au? versus A U9 curve closely follows the theoretical curve up to values of 
A L i 9 = 0.4. Above this value, the restriction that A L37 < 1.2 becomes important, and the curve 
asymptotically tends to the 1 .2 value. For the high absorption bins, Al.t 7 is a constant 1 .2, and 
hence standard deviation envelope goes to zero. 



Fig. 3. The decrease in the observed 37-to-19 GHz absorption 
ratio due to the beamfilling effect. 
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Figure 4 shows the effect of the normalized rms variation (5 on the computation of Al 37 ar >d 
A l)9 . For this figure, equation (35) is used to compute A U 7 and A U9 from the retrieved values 
Al 37 and A L19 using four different p values: 0.7, 0.8, 0.9, and 1.0. That is to say, rather than 
computing (5 for each observation, we use an average value. The theoretical Mie curve lies 
between the P = 0.8 and p = 0.9 curves. This indicates that, on the average, the beamfilling 
effect is characterized by a normalized rms variation P « 0.85, which is somewhat less than the P 
= 1 value given by an exponential probability density function for A L . 



Fig. 4. Parametric curves showing the effect of the normalized 
spatial variability P of liquid water on the 37-to-19 GHz 
absorption ratio. 


8. Inferring Rain Rate from Liquid Water Attenuation 

In this section we describe how rain rate is inferred from the retrieved values of Ali 9 and Al 37 
once the beamfilling correction has been applied. The retrieval method discussed above results 
in A l being directly related to the transmittance x of the radiation from the sea-surface up through 
the atmosphere. Thus, more generally speaking, the retrieved A L is an attenuation factor for the 
polarized sea surface signal traveling from the sea-surface upward through the atmosphere. In 
the absence of scattering, the attenuation and absorption are equivalent and are given by 

H 

A L = | K(h)dh (38) 

0 

where the integral is over altitude h (km), K(h) (km -1 ) is the Mie absorption coefficient, and H 
(km) is the height of the liquid water column. Radiative scattering from rain drops and ice 
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modify the attenuation. For example, the attenuation for point-to-point microwave communication 
links is still given by (38), but K(h) is the Mie extinction coefficient, rather than absorption 
coefficient. However, the sea surface is an extended source, as opposed to a point source. For an 
extended source, the polarized surface signal traveling along the SSM/I viewing direction is 
scattered in other directions, while the surface signals traveling in other directions are scattered 
into the SSM/I viewing direction. These two effects tend to compensate, and using the extinction 
coefficient in (38) would overestimate the attenuation of the surface signal. 

For moderate to heavy rain (R > 10 mm/h), the 19 GHz (37 GHz) extinction coefficient is 
about 20% (60%) higher than the absorption coefficient. One distinguishing characteristic 
between the extinction and absorption coefficients is their spectral signature. For light to 
moderate rain (5 mm/h) the 37-to-19 GHz ratio for the extinction coefficient is 3.8 as compared 
to 3.0 for the absorption coefficient. Figure 3 shows the spectral signature of the SSM/I retrieved 
A L 3 7 /Ali 9 is about 2 for light to moderate rain. Thus, a significantly larger beamfilling correction 
would be needed for the extinction coefficients as compared to the absorption coefficients. We 
decided to use the absorption coefficient to evaluate (38) because 1) its spectral signature is 
closer to the observed Al 37 /A U9i and 2) we expect that the attenuation of the polarized surface 
signal due to scattering will be small (i.e., the scattering in to and out of the viewing direction 
will tend to cancel). 

Fortunately, the choice of the attenuation coefficient does not have a large effect on the 
retrieved rain rate. The larger extinction coefficients would give a lower rain rate except for the 
fact that the beamfilling correction is larger for the extinction coefficients. These two factors 
tend to cancel each other, and, in general, the rain retrievals using the absorption coefficients are 
only about 10% higher than using the extinction coefficients. For example, if the best choice for 
K in (38) is halfway between the absorption and the extinction coefficient, then our rain retrievals 
would be biased about 5% high. 

The Mie absorption coefficient for a rain cloud is found by integrating the Mie absorption 
cross sections (Bom and Wolf, 1975) over the drop size distribution for the rain cloud, as 
described by Wilheit et al. (1977). The drop size distribution is partitioned into two parts: the 
non-precipitating cloud water and the falling rain water. 

K = J N c (r)<r(r)dr + |N R (r)a(r)dr (39) 

where r is the drop radius, N c (r) and N R (r) are the drop size distributions for cloud and rain 
water, respectively, and a(r) is the Mie absorption cross section. When r is much smaller than 
the radiation wavelength, which is the case for the cloud water integral, a(r) is simply 
proportional to r 3 , and as a result the cloud integral is proportional to the total volume of cloud 
water per unit volume of atmosphere. Thus the vertically integrated cloud absorption given by 
(38) is proportional to the columnar cloud water L. 

For the rain integral, the simple a(r) r 3 does not hold, and the absorption depends on the 
details of the drop size distribution. We use the Marshall and Palmer (1948) drop size 
distribution to specify N R (r). The Marshall-Palmer distribution is parameterized in terms of a 
nominal rain rate. Following the method described by Wilheit et al. (1977), we vary this nominal 
rain rate from 0.01 to 50 mm/h and compute the above rain absorption integral, denote by k r , and 
the actual rain rate R assuming the fall velocity given by Waldteufel (1973). We find that the k r 
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versus R relationship in the 19-37 GHz band is well approximated by a power law relationship, 
which is close to linear. 

Combining (38) and (39) gives 

A l ,9 = 0.059(1 - 0.026AT)L + 0.0122(1 + 0.004AT)HR‘ 06 (40a) 

A l37 = 0.208( 1 - 0.026AT)L + 0.0436( 1 - 0.002AT)HR° 95 (40b) 

AT = Tl - 283 (41) 

where L is the columnar cloud water (mm), R (mm/h) is the rain rate averaged over the rain 
column, and H is the height of the rain column. There is a small dependence on the rain cloud 
temperature T L . The mean temperature between the surface and the freezing level is used to 
specify T L . 

T l = (T s + 273)/2 (42) 

where T s is the climatology sea-surface temperature. 

Equation (40) reveals a fundamental problem in retrieving rain rate. Given just Aug and A 1 . 37 . 
it is not possible to uniquely separate and retrieve L, R, and H. The spectral dependencies of the 
cloud water term and the rain rate term are nearly the same, as can be seen by the spectral ratio of 
the coefficients (0.208/0.059 = 3.5;.436/0.0122 = 3.6). And, by doubling the rain rate R and 
halving the height H, one obtains about the same A L . Thus to obtain an estimate for R, one must 
make a priori assumptions regarding L and H. Potentially, these assumptions can produce 
significant errors in the rain retrievals. 

We know little of the statistics of cloud water versus rain rate in rain clouds. The data 
presented by Pruppacher and Klett (1980) indicate that heavy rain clouds can have values of L on 
the order of 1 to 2 mm. It seems reasonable to assume that, in general, L increases with R and 
then tend to level off at very high rain rates. Furthermore, an investigation of 38 northeast 
Pacific extratropical cyclones (Wentz, 1990) indicates that when the SSM/I retrieval of L reaches 
a value of 0.18 mm, a drizzle or light rain is likely. The following expression incorporates the 
features of 1) rain beginning at L = 0.18 mm, 2) L increasing with R, and 3) the L versus R 
relationship leveling off at high R, with L reaching a maximum value between 1 and 2 mm. 

L = 0.18 (l + VHR) (43) 

The specification of the rain column height H is based, in part, on the altitude H F of the 
freezing level as derived from the radiosonde observations. The global radiosonde observations 
for the 1987-1990 period collected by Wentz (1997) are used to find H F as a function of the sea- 
surface temperature T s . Out of the total 42,195 radiosonde soundings, we only use the 9,120 
soundings for which the surface relative humidity is > 90%. By restricting the data set to high 
humidity cases, the results should be more indicative of rain observations. Figure 5 shows the 
height of the freezing level measured by the radiosondes versus the climatological sea surface 
temperature at the radiosonde site. For the stations at very high latitudes, the typical value of H F 
is about 1 km. The mid-latitude value of H F ranges from 2 to 4 km, and in the tropics H F reaches 
a value of 5 km. 

Equation (40) shows that the retrieved rain rate is very nearly proportional to H _l . Thus the 
proper specification of H is critical to obtaining good rain rate retrievals. In a preliminary 


18 



analysis, we used the H F values shown in Figure 5 to specify H and found that the rain rates in 
the tropics were about 40% low in comparison to other climatologies (see Section 9). We 
attribute at least part of this systematic underestimation to cases for which the rain column does 
not extend the full 5 km to the freezing level. In general, H F is an upper bound on H, and there 
will be many cases for which H < H P . Thus using H F will, on the average, tend to underestimate 
the rain rate. If there is a substantial amount of rain in the tropics that forms without an ice phase 
(Fletcher, 1969) then the underestimation will be large. 

In the tropics, we find that using H = 3 km, rather than the 5-km freezing level, corrects the 
underestimation of rain relative to other climatologies. We believe that the assumption of a mean 
column height of 3 km is not unreasonable and may be an indication of the prevalence of warm 
rain processes in the tropics. Outside the tropics, we use H F to specify H. The following simple 
regression, which is shown in Figure 5, is derived so as to match H F for low values of Ts and to 
reach a value of 3 km for high values of Ts. 

H = 1 + 0. 1 4(T S - 273) - 0.0025(T S - 273 f T s < 30 1 (44a) 

H = 3 T s > 301 (44b) 



Fig. 5. The altitude of the atmospheric freezing level plotted 
versus the climatology sea surface temperature. 
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Having specified H and the relationship between R and L, one can invert equation (40) and 
find a value for R given A L . Note that as a result of the beamfilling correction discussed in 
Section 7, the retrieved values of A L i 9 and A L37 are not independent. Rather they are computed 
such that their ratio is consistent with equation (40) above. For this reason, the same value for R 
is found from either Ali 9 or Al 3 7 . The one exception is when Al 3 7 exceeds the maximum value 
of 1.2. In this case, A L | 9 is used to compute the rain rate. It is important to note that retrieved 
value of R is the vertically averaged rain rate given by 

R = H _1 fdh R(h) (45) 

Jo 

where R(h) is the rain profile. The difference between R and R(0) is an additional source of error 
when comparing to in situ surface rain measurements. Evaporation under the rain cloud will tend 
to make R greater than R(0), while a decrease in R(h) at the top of the rain cloud will tend to 
make R less than R(0). 

9. Rain Retrieval Results 

а. Probability density function of SSM/1 rain rates 

All results in this section are based on SSM/I observations for the four year period from 1991 
through 1994. Observations from two satellites, F10 and FI 1, are used. The F10 observations 
cover the entire 4 year period, while the FI 1 observations begin in January 1992. The top frame 
in Figure 6 shows the probability density function (pdf) for the rain rates inferred from the two 
SSM/I’s. The thick curve shows global results, and the thin curve shows tropical results (20S to 
20N). The computation of any rain pdf is very dependent on the temporal and spatial averaging. 
For the SSM/I, the temporal averaging is essentially instantaneous, and the spatial averaging has 
a resolution of about 32 km. A rain pdf computed from rain gauges looks very different than that 
shown in Figure 6 because the spatial averaging is very different. The leftmost point on the pdf 
curves corresponds to the number of no-rain observations. A total of 85.9% of the SSM/I 
observations indicated no rain. An additional 8.3% of the observations indicated very light rain 
not exceeding 0.2 mm/hr, and the remaining 5.8% of the observations indicate rain exceeding 
0.2 mm/hr. We consider the accuracy of the ‘very light rain’ retrievals as questionable. Some or 
many of these observations may actually be heavy non-raining clouds. Note that the contribution 
of the very light rain observations to the total rainfall is very small (see below). 

To determine the contribution of the various footprint-averaged rain rates to the overall 
rainfall amount, we multiply the rain pdf by the rain rate, as shown in the bottom frame in Figure 

б. In this case, the area under the curve equals the average oceanic rainfall, which is 0.12 mm/hr 
(2.9 mm/day) globally and 0.16 mm/day (3.9 mm/day) in the tropics. The questionable very light 
rain observations (R < 0.2 mm/h) only contribute 0.007 mm/h (0.17 mm/day) to this total. One- 
half of the total global oceanic rainfall occurs at footprint-averaged rates above (and below) 
about 3.5 mm/h. For rainfall in the tropics, this value increases to 5.5 mm/h. Due to the large 
size of the footprint (32 km) over which the enveloped rainfall is averaged, this mid-point value 
is much lower than that obtained from rain gauges. Four-minute rain gauge statistics (Jones and 
Sims, 1 978) suggest that about half of tropical rainfall occurs at rates above about 20 mm/h. One 
possible interpretation of this result is that, on the average when significant rain is being 
observed, only about one quarter of the SSM/I footprint is actually covered by rain. 
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Fig. 6. The top frame show probability density function (pdf) for rain rate averaged 
over the SSM/I footprint. The bottom frame show the rain pdf times the rain rate. 
The thick curves show global results, and the thin curves show tropical results. 


b. Global distribution of SSM/I rain rates 

Figure 7 shows the seasonal and annual zonally averaged rainfall computed from the SSM/I 
observations for 1991 through 1994. The meridional structures revealed by the SSM/I are similar 
to previously published climatologies. The maximum oceanic rainfall occurs at the equatorial 
latitudes associated with the strong convection in the intertropical convergence zone (ITCZ) for 
all seasons. This peak is quite narrow in meridional extent and varies from about 7 mm/day in 
the winter to a maximum 1 1 mm/day in the summer. The seasonal north-south migration of the 
ITCZ, which is in phase with the solar insolation, is also apparent in the figure. The extratropical 
rainfall is greater in the northern hemisphere than in the southern hemisphere for all seasons. 
Low precipitation rates (~ 1 mm/day) are observed in those zones of subsidence influenced by 
the large semipermanent anticylones. 

Figure 8 shows color-coded global maps of the SSM/I annual and seasonal rainfall average 
over the four years from 1991 through 1994. The major features of the spatial distribution of the 
average annual rainfall are quite similar to those revealed in other satellite climatologies (see 
below). The largest annual rainfall amounts are seen to occur in the tropical Pacific, extending 
from South America to Papua New Guinea. Peaks of 15 mm/day occur throughout this band. 
Additional heavy rain associated with the Indian summer monsoons is apparent in the Bay of 
Bengal. The other major feature of the global rainfall maps is the extremely dry areas associated 
with the large semipermanent anticylones in the Southeast Pacific and Southeast Atlantic. These 
areas are essentially void of rain (R < 0.3 mm/day). 
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Fig. 7. The zonal average of rainfall over the ocean. The thick curve, which is 
repeated in each frame, is the annual average. The thin curve is the 3-month 
seasonal average. 


c. Comparison to other satellite climatologies 

We now compare our rainfall estimates (hereafter “WS”) to two other emission-based rain 
climatologies: Spencer (1993; hereafter “MSU”), and Wilheit et al. (1991; hereafter “WCC”). 
The MSU rain rates are inferred from the 50.3 GHz T B observations taken by the Microwave 
Sounding Unit (MSU). The WCC rain rates are inferred from the SSM/I T B observations. The 
same period of record (1991-1994) is used from these data sets. Figure 9 compares the three 
estimates of the annual zonally averaged rainfall. In general, the three rainfall estimates are 
similar, but there are some notable differences. We first note that above 50N and below 55S, the 
MSU rain data are contaminated by sea ice (see below). This explains the upturn at the two ends 
of the MSU curve in Figure 9. In the ITCZ, the WS, MSU, and WCC reach maximum values of 
8.1, 7.4, and 6.9 mm/day, respectively. This represents about a 15% difference between the 
highest estimate (WS) and the lowest estimate (WCC). In the extratropics storm track regions, 
the situation changes. Here the WS rainfall is the lowest, and MSU is the highest. Very close 
agreement is seen in the very dry areas associated with the semipermanent anticylones. 

Figure 10 shows color-coded global maps of the MSU minus WS rainfall and the WCC minus 
WS rainfall. To compute these differences, the rainfall is averaged over the four years (1991- 
1994) and then smoothed to a spatial resolution of about 300 km. The largest differences are 
seen between the MSU and WS. The MSU produces more rainfall in the downstream portions of 
the extratropical storm tracks, and less rainfall over most portions of the tropics, particularly in 
the tropical west Pacific. Comparisons of Figure 10 to SSM/I retrievals of cloud water (not 
shown) suggest that the MSU-WS differences might be related to cloud water. Areas where the 
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Annually and Seasonally Averaged Rainfall (mrrVday) 


Fig. 8. Average rainfall for 1991 through 1994 derived from SSM/I. The top 
map shows the annual rainfall, and the four smaller maps show the seasonal 
averages. 

MSU-WS difference is significantly positive (negative) are moderately correlated with areas 
having a relatively high (low) cloud amount as compared to the rainfall. One example Ls the 
downstream portions of the extratropical storm tracks where there is significant cloud coverage 
but relatively little rain. In these regions the MSU rainfall is about 2 mm/day higher than WS. In 
contrast, along most of the ITCZ, the cloud content is relative small compared to the heavy rain, 
the MSU rainfall is about 2 to 3 mm/day lower than WS. An interesting ocean area is seen just 
west of Central America and Columbia. The north (south) part of this area shows large negative 
(positive) MSU-WS differences. An analysis of SSM/I retrievals shows moderately heavy rain 
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Fig. 9. A comparison of zonally averaged rain rates from three satellite 
climatologies. WS denotes our results, WCC denotes the Wilheit et al. 
(1991) results, and MSU denotes the Spencer (1993) results. The upturn 
at the two ends of the MSU curve is due to sea ice contamination. 


and relatively small cloud contents in the north and the reverse situation in the south, which is the 
same correlation as seen in the storm tracks and the UCZ. The correct partitioning of cloud and 
rain water is a problem for both MSU and SSM/I. As pointed out by Spencer (1993), the 
hypersensitivity of the MSU 50.3 GHz channel to both cloud water and rain water makes the 
MSU unable to distinguish between the two. We have somewhat more confidence in the SSM/I 
rainfall because the frequencies of 19.3 and 37 GHz are less sensitive to cloud water, and we 
have attempted to do a cloud versus rain partitioning. This confidence is bolster by the fact that 
the cloud-to-rain ratio derived from SSM/I seems realistic. It is a minimum just off the east 
coasts of the continents where baroclinic wave activity is the strongest. Then this ratio increases 
eastward across the ocean basins, consistent with weaker wave activity. 


The difference map between WS and WCC shows better agreement. The major difference is 
in tropical areas of heavy rain, where the WS is about 2 mm/day higher. In the extratropical 
storm tracks, the WCC is typically about 1 mm/day higher. In the dry areas, all three rain 
estimates (WS, MSU, and WCC) agree well. We find no obvious correlation between the WS- 
WCC difference and other parameters, except for the rainfall itself. When the rain is very heavy, 
WS tends to be higher than WCC. 

Note that in the MSU-WS figure, the red areas in Sea of Okhotsk, the Bering Sea, Hudson 
Bay, Labrador Sea, and off Antarctica are sea ice contamination in the MSU rain product. A very 
small amount of ice contamination is also seen in the WCC product just north of Japan. 
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Difference in Rainfall for the 1991 Through 1994 Period (mnVday) 


Fig. 10. A comparison of three satellite-derived rain maps. The top image shows the 
rainfall derived from the MSU (Spencer 1993) minus the SSM/I rainfall computed 
from the algorithm described herein (WS). The bottom image shows the rainfall 
produced by the Wilheitet al. (1991) algorithm (WCC) minus the WS rainfall. 


10. Conclusions 

A new method for the physical retrieval of rain rates from the SSM/I has been presented and 
compared to two other microwave emission-type rainfall climatologies. The method is part of a 
unified ocean parameter retrieval algorithm that also diagnoses total integrated water vapor, cloud 
water, and wind speed. The algorithm first retrieves those geophysical parameters directly 
observed by SSM/I. With respect to rain, the two relevant parameters are the liquid water 
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transmittance x L and its corresponding effective radiating temperature Ty. A beamfilling 
correction is then applied based on the spectral signature of the retrieved x L at 19 and at 37 GHz. 
Finally assumptions about the cloud versus rain water partitioning and the rain column height are 
made to obtain a value for rain rate. 

We find that the water vapor retrievals maintain reasonably good accuracy when there is rain 
in the field of view. The rms difference between the SSM/I water vapor retrieval and 
radiosondes is about 5 mm for rain rates from 1 to 15 mm/h. Errors in the water vapor retrievals 
appear to be uncorrelated with the rain rate. This indicates that the algorithm is successfully 
orthogonalizing the retrievals so that there is a minimum of cross-talk between the retrieved 
parameters. 

As expected the effective radiating temperature Ty exhibits a strong depression relative to the 
mean air temperature for moderate to heavy rain. This depression is due 1) to radiative scattering 
from large rain drops and ice and 2) to that fact that most of the radiation is coming from the cold 
top part of the rain cloud. For the heaviest rain, the Ty depression is -10 K and -20 K for 19 and 
37 GHz, respectively. 

The spectral signature of the retrieved x L shows that the ratio of the 37-GHz to 19-GHz liquid 
water absorption is, on the average, about 40% lower than predicted by Mie theory for moderate 
to heavy rain. We attribute this difference to the beamfilling effect, which we parameterize in 
terms of the normalized rms variation (3 of the liquid water absorption A L . To correct for this 
effect, the 37-GHz to 19-GHz liquid water absorptions are increased until the Mie ratio is 
realized. Globally, we find the normalized rms variation of Al to be about 0.85, which is 
somewhat less than that for an exponential pdf. Of course, individual observations can have 
values of P substantially higher than 0.85. 

The last step in the rain rate retrieval involves dividing the columnar rain rate by some 
assumed rain column height H. In the tropics, we find using the freezing level, which is about 
5 km, to specify H results in tropical rain rates that appear to be to low when compared with 
other rainfall climatologies. This low bias may be due to a substantial portion of the rain in the 
tropics not extending all the way to the freezing level (Fletcher, 1969). To correct the low bias, 
we use a value of H ~ 3 km in the tropics. 

Global rain rates are produced for the 1991 through 1994 period from two SSM/I’s that are 
aboard the F10 and FI 1 satellites. The probability density function (pdf) for the retrieved rain 
rates (averaged over the 32 km SSM/I footprint) shows that on a global basis 6% of the SSM/I 
observations detect measurable rain rates (R > 0.2 mm/h). One-half of the total global oceanic 
rainfall occurs at footprint-averaged rates above (and below) about 3.5 mm/h. For rainfall in the 
tropics, this value increases to 5.5 mm/h. Globally, the average rainfall over the oceans is about 
2.9 mm/day, and in the tropics (20 N to 20 S) it is 3.9 mm/day. Zonal averages and global maps 
of the retrieved rain rates show structures that are similar to those in previously published rain 
climatologies. The peak at the ITCZ is quite narrow in meridional extent and varies from about 
7 mm/day in the winter to a maximum 1 1 mm/day in the summer. Very low precipitation rates 
(< 0.3 mm/day) are observed in those areas of subsidence influenced by the large semipermanent 
anticylones. 

We compare our rainfall estimates (WS) to two other emission-based rain satellite 
climatologies: Spencer (1993; MSU) and Wilheit et al. (1991; WCC). The agreement is 


26 



reasonably good. In the ITCZ, the WS, MSU, and WCC reach maximum values of 8.1, 7.4, and 
6.9 mm/day, respectively. In the extratropics storm track regions, the situation reverses, with the 
WS rainfall as the lowest, and MSU as the highest. Very close agreement is seen in the very dry 
areas associated with the semipermanent anticylones. The largest differences are seen between 
the MSU and WS. In the downstream portions of the extratropical storm tracks, the MSU rain is 
significantly greater than the WS. In contrast, the MSU rainfall is less over most portions of the 
tropics, particularly in the tropical west Pacific. The MSU minus WS rain difference seems to be 
related with cloud water. Areas where the MSU-WS difference is significantly positive 
(negative) are moderately correlated with areas having a relatively high (low) cloud amount as 
compared to the rainfall. 

Our rain retrieval technique could probably be improved by including the SSM/I 85 GHz 
channels. These channels are very sensitivity to radiative scattering by ice and may provide the 
means to better identify areas of heavy rain exceeding 15 mm/h. 

In conclusion, there still remains the problem of absolutely calibrating the rain algorithm. The 
lack of good-quality in situ rain measurements over the oceans has been a major source of 
difficulty for all satellite -based rainfall estimation techniques, and it is still not clear how to best 
deal with the calibration problem. Hopefully future programs such as TRMM and the 
Precipitation Intercomparison Project (PIP) will contribute to the better calibration of rainfall 
derived from satellites. 


Acknowledgments. This research was supported by NASA’s Oceans Program and EOS 
Program under contracts NASW-4714 and NAS5-32594. We are thankful to the Defense 
Meteorological Satellite Program for making the SSM/I data available to the civilian community. 


27 



REFERENCES 

Buettner, K.J.K., 1963: Regenortung vom wettersatelliten mit hilfe von zentimeterwellen (Rain 
localization from a weather satellite via centermeter waves). Naturwiss., 50 , 591 . 

Fletcher, N.H., 1969: The Physics of Rainclouds. Cambridge University Press, New York, NY, 
390pp. 

Goldstein, H., 1951: Attenuation by condensed water. Propagation of Short Radio Waves, MIT 
Rad. Lab. Ser., 13, McGraw-Hill, New York. 

Hollinger, J., R. Lo, G. Poe, R. Savage, and J. Pierce, 1987: Special Sensor Microwave/Imager 
User’s Guide. NRL Tech. Rpt., Naval Research Laboratory, Washington, DC. 120 pp. 

Jones, D.M.A., and A.L. Sims, 1978: Climatology of instantaneous rainfall rates. J. Appl. 
Meteor., 17 , 1 135-1 140. 

Petty, G. W., 1994: Physical retrievals of over-ocean rain rate from multichannel microwave 
imagery. Part I: Theoretical characteristics of normalized polarization and scattering 

indices. Meteorol. Atmos. Phys., 54 , 79-99. 

Shea, D.J., K.E. Trenberth, and R.W. Reynolds, 1990: A global monthly sea surface temperature 
climatology. NCAR Tech. Note 345, National Center for Atmospheric Research, 167 pp. 

Simpson, J., R.F. Adler, and G.R. North, 1988: A proposed Tropical Rainfall Measuring Mission 
(TRMM) satellite. Bull. Amer. Meteor. Soc., 69 , 278-295. 

Spencer, R.W., 1986: A satellite passive 37 GHz scattering-based method for measuring oceanic 
rain rates. J. Climate Appl. Meteor., 25, 754-766. 

Spencer, R.W., 1993: Global oceanic precipitation from the MSU during 1979-92 and 
comparisons to other climatologies. J. Climate, 6, 1301-1326. 

Spencer, R.W., H.M. Goodman, and R.E. Hood, 1989: Precipitation retrieval over land and 
ocean with the SSM/I: identification and characteristics of the scattering signal. J. Atmos. 
Oceanic Tech., 6, 254-273. 

Wentz, F.J., 1990: SBIR Phase II Report: West coast storm forecasting with SSM/I, RSS Tech. 
Rpt. 033190, Remote Sensing Systems, Santa Rosa, Ca, 378 pp. 

Wentz, F.J., 1991: User’s manual: SSM/I antenna temperature tapes, revision 1, RSS Tech. Rpt. 
120191, Remote Sensing Systems, Santa Rosa, CA, 69 pp. 

Wentz, F.J., 1997: A well-calibrated ocean algorithm for SSM/I, J. Geophy. Res., in press. 

Wilheit, T.T., A.T.C. Chang, M.S.V. Rao, E.B. Rodgers, and J.S. Theon, 1977: A satellite 
technique for quantitatively mapping rainfall over the oceans. J. Appl. Meteor., 16 , 551-560. 

Wilheit, T.T., A.T.C. Chamg, and L.S. Chiu, 1991: Retrieval of monthly rainfall indices from 
microwave radiometric measurements using probability distribution functions. J. Atmos. 
Oceanic Tech., 8, 1 18-136. 


28 



NASA 

Nnioui Acavmics and 


1. Report No. 


4. Title and Subtitle 

SSM/I Rain Retrievals 
Ocean Algorithm 


Report Documentation Page 


2 . 


Government Accession No. 


3. Recipient's Catalog No. 


5. Report Date 


Within A Unified All-Weather 


December 1, 1996 

6. Performing Organization Code 


7. Author(s) 

Frank J. Wentz 
Roy W. Spencer 


9. Performing Organization Name and Address 
Remote Sensing Systems 
1101 College Avenue, Suite 220 
Santa Rosa, CA 95404 


12. Sponsoring Agency Name and Address 

NASA 

Mission to Planet Earth 
NASA HQ, Code YS 
Washington, DC 20546 

15. Supplementary Notes 


8. Performing Organization Report No. 

RSS Tech, Report 120196 

10. Work Unit No. 


11. Contract or Grant Np. 

NASW - 4714 

13. Type of Report and Period Covered 

Final Report, 1992-1996 

14. Sponsoring Agency Code 


A new method for the physical retrieval of rain rates from satellite microwave radiometers is presented and compared to two 

other rainfall climatologies derived from satellites. The method is part of a unified ocean parameter retrieval algorithm that is based 

16. Abstract 0 n the fundamental principles of radiative transfer. The algorithm simultaneously finds near-surface wind speed W, columnar water 
vapor V, columnar cloud liquid water L, rain rate R, and effective radiating temperature Tu for the upwelling radiation. The 
performance of the algorithm in the absence of rain is discussed in Wentz [1997], and this paper focuses on the rain component of the 
algorithm. A particular strength of the unified algorithm is its ability to ‘orthogonalize’ the retrievals so that there is minimum cross- 
talk between the retrieved parameters. For example, comparisons of the retrieved water vapor with radiosonde observations show that 
there is very little correlation between the water vapor retrieval error and rain rate. For rain rates from 1 to 15 mm/h, the rms 
difference between the retrieved water vapor and the radiosonde value is 5 ram. A novel feature of the rain retrieval method is a 
beamfilling correction that is based upon the ratio of the retrieved liquid water absorption coefficients at 37 GHz and 19.35 GHz. 

This ratio decreases by about 40% when heavy and light rain co-exist within the SSM/I footprint as compared to the case of uniform 
rain. This correction has the effect of increasing the rain rate when the spectral ratio of the absorption coefficients is small. Even 
with this beamfilling correction, tropical rainfall is still unrealistically low when the freezing level in the tropics (~ 5 km) is used to 
specify the rain layer thickness. We restore realism by reducing the assumed averaged tropical rain layer thickness to 3 km, thereby 
accounting for the existence of warm rain processes in which the rain layer does not extend to the freezing level. Global rain rates are 
produced for the 1991 through 1994 period from observations taken by microwave radiometers (SSM/I) that are aboard two polar- 
orbiting satellites. We find that approximately 6% of the SSM/I observations detect measurable rain rates (R > 0.2 mm/h). Zonal 
averages of the rain rates show the peak at the intertropical convergence zone (ITCZ) is quite narrow in meridional extent and varies 
from about 7 mm/day in the winter to a maximum 11 mm/day in the summer. Very low precipitation rates (< 0.3 mm/day) are 
observed in those areas of subsidence influenced by the large semipermanent anticylones. In general, these features are similar to 
those reported in previously published rain climatologies (Wilheit et al., 1991; Spencer, 1993). However, significant differences do 
exists between our rain rates and those produced by Spencer (1993). These differences seem to be related to non-precipitating cloud 
water. 


17. Key Words (Suggested by Author(s)) 
Microwave radiometer 
Precipitation 


18. Distribution Statement 


19. Security Classif. (of this report) 


20. Security Classif. (of this page) 


21. No. of pages 

28 


22. Price 


NASA FORM 1626 OCT 86 




















